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Abstract 

Recent observation of pulsar PSR J1614-2230 with mass about 2 solar masses poses a severe 
constraint on the equations of state (EOS) of matter describing stars under extreme conditions. 
Neutron stars (NS) can reach the mass limits set by PSR J1614-2230. But stars having hyperons or 
quark stars (QS) having boson condensates, with softer EOS can barely reach such limits and are 
ruled out. QS with pure strange matter also cannot have such high mass unless the effect of strong 
coupling constant or color superconductivity are considered. In this work I try to calculate the 
upper mass limit for a hybrid stars (HS) having a quark-hadron mixed phase. The hadronic matter 
(having hyperons) EOS is described by relativistic mean field theory and for the quark phase I 
use the simple MIT bag model. I construct the intermediate mixed phase using Glendenning 
construction. HS with a mixed phase cannot reach the mass limit set by PSR J1614-2230 unless 
I assume a density dependent bag constant. For such case the mixed phase region is small. The 
maximum mass of a mixed hybrid star obtained with such mixed phase region is 2.01Mq. 
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I. INTRODUCTION 



Neutron stars (NS) are gravitationally bound, therefore the precise measurement of mass 
and radius of a NS should provide a very fine probe for the equation of state (EOS) of dense 
matter. The first reasonable ideas about the composition of compact stars argued that 
matter are under extreme densities and is mainly composed of neutrons with small fractions 
of protons and electrons. Further theoretical developments and modern experimental results 
opened the window to other possibilities. The densities in the interior of neutron stars is 
about 3 — 10 times that of the nuclear saturation density {no ~ 0.15 fm~^). At such high 
densities in their interiors, the matter there is likely to be in a deconfined and chirally 
restored quark phase |l|. 

The strange matter hypothesis was first proposed by Itoh and Bodemer and was 

then improved by Witten |^. It states that, matter at extreme density and/or temperature 
are composed of almost equal number of up, down and strange quarks, called strange quark 
matter (SQM). It is also the ground state of strongly interacting matter at such extreme 
conditions. If this is true, then matter at such extreme conditions is likely to eventually 
convert to SQM. Such a high density scenario is present in the interiors of a NS and therefore 
normal nuclear matter is likely to undergo a phase transition and converts to SQM. The 
strange matter hypothesis was first extensively studied in the simple MIT bag model by 
Farhi & Jaffe The conversion process and the phase transition was further analyzed by 
Alcock et al. [6]. The phase transition in a NS may continue up to the surface of the star 
or may stop inside the star. Depending up on where this, a quark star (QS) may be of 
two types, a strange star (SS) or a hybrid star (HS). SS are stars composed only of SQM, 
while HS has a quark core and a hadronic exteriors. IN the region between the quark core 
and hadronic outer matter, there may exist a mixed phase region where both quarks and 
hadrons are present. Thus the observed pulsars are still very much model dependent. 

Recently, Demorest et al. [7] found a new maximum mass limit for compact stars by 
measuring very precisely the mass of the millisecond pulsar PSR J16 14-2230 to be M = 
1.97 ± 0.04 Mq. This value is much higher than any previously measured pulsar mass. This 
measurement, has imposed a very severe constraints on the EOS of matter describing the 
compact objects. The model of NS, without hyperons, can easily satisfy the new mass 
constraint. However the presence of strangeness, either in the form of hyperons in nuclear 



matter or in the form od strange quarks in quark matter, cannot easily satisfy the mass 
hmit. So, new studies had been carried out to make the hyperonic EOS and quark EOS to 
satisfactorily explain the new mass constraint. 

Basically to satisfy the new mass limit, one has to make the EOS stiffer, which usually 
is softened by the presence of strangeness. In the hyperonic nuclear matter sector, recent 
studies have suggested that the stiffening of hyperonic EOS is possible at par with the 
new experimental results 8|. Authors also had revisited the role of vector meson-hyperon 
coupling {9 1 and hyperon potentials 10|, to calculate the maximum mass. 

Studies prior to the discovery pulsar PSR J1614-2230 have suggested the stiffening quark 
matter EOS from the effect of strong interactions, such as one-gluon exchange or color- 
superconductivity llNl7l|. which can satisfy the new constraint. Ozel 18| and Lattimer 



19| gave first studies on the implications of the new mass limits from PSR J1614-2230 for 



quark and hybrid stars in the quark bag model. Recently, Bonanno & Sedrakian [20[ has 
succeeded in obtaining massive HS. They employed color-superconducting quark core and 
very stiff hadronic EOS (like the NL3 hyperonic model or the GM3 nuclear model). 

In this work I perform an extensive study of hybrid star mass using the relativistic mean- 
field hadronic EOS together with a simple three-flavor MIT bag model quark EOS. The 
model of the HS has a mixed phase intermediate region. I would also discuss as how the 
understanding of more precise astrophysical measurements of the mass and radius of neutron 
stars can help revealing the viability of exotic quark star models. The paper is organized 
as follows: In Section II, I describe the hadronic phase and in section III I describe the 
MIT bag model. The mixed phase EOS is constructed in section IV, using the Glendenning 
construction. I present my plots and extensively describe my results for the EOS and the 
mass-radius curve in Section V. The maximum mass for the hybrid star is also calculated in 
this section. Finally in section VI, I summarize my results and draw important conclusion 
from them. 



HADRONIC PHASE 

At the outermost region of the star, at comparatively low densities the matter is mainly 
composed of hadrons. I use the non linear relativistic mean field (RMF) model with hyperons 
(TMl parametrization) to describe the hadronic phase EOS. In this model the baryons 
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interact with mean meson fields 2l|-|25|. 

The model lagrangian density includes nucleons, baryon octet (A, T,^'"^, and leptons 

b ^ 

- (mb - gabiy)]^pb + -{df,ad^(T - mla^) 

11 11 



(1) 



Leptons L are non-interacting but the baryons are coupled with the scalar a mesons, the 
isoscalar-vector mesons and the isovector-vector mesons. The model constants are 
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25|. The 



fitted according to the experimental results of bulk properties of nuclear matter 
TMl model explains the nuclear saturation of but cannot sufficiently models the hyperonic 
matter, as it fails to reproduce the strong observed AA attraction. This defect can be 
remedied by Mishustin & Schaffner 25| by the addision of iso-scalar scalar a* mesons and 
the iso-vector vector mesons, coupling only with the hyperons. 

The detailed EOS calculation can be found in the above mentioned references 
and I do not repeat them here. 

The total energy density takes the form 
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and the pressure can be represented as 



(2) 



P = ^Pini-e, (3) 

i 

where /ij and is the chemical potential and number density of particle species i. 



QUARK PHASE 



The quark phase is modeled according to the simple MIT bag model |26|. The cur- 
rent masses of up and down quarks are extremely small, e.g., 5 and 10 MeV respectively, 



whereas, for strange quark the current quark mass is not well established, and I vary it in 
my calculation. For the bag model the energy density and pressure can be written as 

E ^ f'dkk'^m^ + k^ + Ba, (4) 

i=u,d,s ^ 

E ^2 dk ' -Bo, (5) 

where kp = ^Jfij — mf and gi is the Fermi momentum and degeneracy factor of quarks of 
species i. Bq is the energy density difference between the perturbative vacuum and the true 
vacuum, i.e., the bag constant. In this sense Bq can be considered as a free parameter. 

Both the hadronic and quark matter, maintains baryon number conservation, and are 
beta-equilibrated and charge neutral. 




MIXED PHASE 



With the previously described hadronic and quark EOS, Glendenning construction [27 1 
gives the mixed phase regime. The mixed phase is the baryong density range where both 
quarks and hadrons are present. In the mixed phase the hadron and the quark phases are 
separately charged but the mixed phase is charge neutral as a whole. Thus the matter can 
be parametrized by the pair of electron and baryon chemical potentials /^e and yU„. Pressure 
of the two phases are made equal to maintaining mechanical equilibrium. To satisfy the 
chemical and beta equilibrium conditions the chemical potential of different particles are 
related to each other. The Gibbs criterion gives the mechanical and chemical equilibrium 
between two phases, and is written as 

Pnv{.l^ej l^n) = -PQp(/^e5 /^n) = -PfvIP- (6) 

The solution of above equation gives the equilibrium chemical potentials of the mixed phase. 
As the two phases intersects one can calculate the corresponding charge densities of the 
hadronic components p^^ and quark components p^^ separately in the mixed phase. The 
volume fraction occupied by quark matter in the mixed phase x is given by 

XpT + (1 - x)pT = 0. (7) 
The mixed phase energy density eMP and the number density tt-mp can be written as 

CMP = Xeqp + (1 - x)eHP, (8) 
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FIG. 1. Pressure as a function of energy density with bag pressure of 170 and ISOMeV. 
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FIG. 2. Pressure as a function of baryon density with bag pressure of 170 and ISOMeV. 



^MP = X^Qp + (1 - x)nHp. 



(9) 



Therefore the EOS is now a system having a charge neutral hadronic phase at lower 
densities, a charge neutral mixed phase in the intermediate region and a charge neutral 
quark phase at higher densities. 

RESULTS 

The EOS are constructed to describe the properties of matter inside a NS, therefore the 
EOS properties would also resemble the properties of a NS. The central region of the star 
has maximum density (few times no) , therefore the matter at the core is most likely to have 
a phase transition. Therefore the central region would have stable strange matter (or a 
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colour superconducting matter). As the density decreases radially outwards some nuclear 
matter (nucleons) starts appearing and so in the intermediate region there is likely to have 
a mixed phase. Much further outwards I have only matter consisting of only nucleons. The 
crust consisting mainly of free electrons and nuclei, which completes the star structure. 

The hadronic EOS, I assume a fixed TMl parameter set, which satisfactorily explains 
the properties of hadronic matter at extreme condition. I can control the quark EOS by 
changing the masses of the quarks and the bag constant. The masses of the light quarks are 
quite bounded and take them to be 5MeV (u) and lOMeV (d). The mass of s-quark is still 
not well established, but expected to lie between 100 — SOOMeV. I would vary the mass of 
the s-quark within this bounded mass range. I would also vary the bag constant (Be) to 
regulate the mixed phase region. This parametrization of the EOS of the hadron and quark 
matter is responsible for characterization of the matter in the mixed phase region. Using 
the Glendenning construction to construct the mixed phase, and plot curves of pressure 
against energy density as seen in fig [H In fig [T] I have plotted the mixed phase EOS with 
bag pressures 170MeV and ISOMeV. Actually the relation runs as Bg^/^ = 170MeV, but for 
simplicity I will denote Bg^^^ = 170MeV = Bg. For this case the mass of the s-quark (m^) 
is taken to be ISOMeV. With a constant bag pressure, lower bag pressure cannot generate a 
mixed phase region. I do not go above Bg = ISOMeV, as for that case the EOS becomes very 
flat, and the maximum mass of the star becomes less. In the curves, the lower portion is 
nuclear phase (dotted/dash line), the intermediate region is the mixed phase (bold hne) and 
the higher region is the quark phase (dotted/dash line). Fig[2l shows the pressure against 
baryon density for bag constant 170MeV and ISOMeV. The mixed phase starts at 0.2/m~^ 
and ends at 0.76/m~^ for bag pressure 170MeV. For bag pressure ISOMeV the mixed phase 
region is in between 0.22fm~^ and O.S9/m~^. The curve with bag constant 170MeV is much 
stiffer than the curve with bag pressure ISOMeV, because the bag pressure adds negatively 
to the matter pressure, making the effective pressure low. The above curves also shows 
that as the bag pressure increases the range of mixed phase region also increases. As the 
variation of pressure with both energy density and baryon density is quite similar, from now 
on I would only plot curve showing pressure as function of energy density. 

With such high bag pressure it is impossible to attain the mass limit set by PSR J1614- 
2230. Therefore I have to devise some other mechanism which would give stiffer EOS, 
thereby increasing the maximum mass of the HS. For that, I assume a density dependent bag 
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FIG. 3. Pressure against energy density plot with constant and varying bag pressure, having 



Bg = ITOMeV. 



constant. In the literature there are several attempts to understand the density dependence 
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29| : however, currently the results are highly model dependent and still there is no 



definite picture. I parametrized the bag constant in such a way that it attains a value -Boo, 
asymptotically at very high densities. The range of value of obtained from experiments 
can be found in Burgio et al. 30| , and I assume it to be 130MeV, the lowest value mentionec 



there. With such assumptions I then construct a Gaussian parametrization given as 



30 



3l| 



Bgniub) = Boo + (Bg - 5oo) exp 



-/3 



no' 



(10) 



The lowest value of Bgn, which is its value at the asymptotic high density in quark matter, 
is fixed at ISOMeV. The bag pressure quoted would be the value of the bag constant at 
the starting of the mixed phase region on the low density regime {Bg in the equation). As 
the density increases the bag pressure decreases and reaches ISOMeV asymptotically, the 
decrease rate is controlled by (3. 

In fig [3] I have plotted curves showing the difference in the slope of the curves with and 
without the variation of bag pressure (for Bg = ITOMeV). For the varying bag pressure 
the mixed phase region shrinks, becomes flatter but the quark phase region becomes stiffer. 
The mixed phase region now only extends up to baryon density 0.53/m~^. The change 
in the mixed phase region is about ~ 30%. This is because, going to higher densities (or 
higher energy density towards the core) the effective matter pressure increases with the 
decrease in bag pressure (bag pressure adds negatively to the matter pressure). With such 
a density dependent bag constant I can have a signiflcant mixed phase region with lower 
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FIG. 4. Pressure against energy density plot with varying bag pressure, having Bg = 160 and 
ISOMeV. 
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FIG. 5. Pressure vs energy density plot showing the explicitly the mixed phase region, for the 
varying bag pressure Bg = 150MeV. 

values of bag pressure. As shown in fig H] I can have mixed phase region with bag pressure 
Bg, for leOMeV and ISOMeV. For the 160MeV EOS the s-quark mass = ISOMeV and 
for the ISOMeV curve the s-quark mass is = 300MeV. With bag pressure, Bg, 160 and 
ISOMeV the mixed phase region is of considerable small. For bag constant 160MeV the 
mixed phase region starts at density 0.1S/m~^ and ends at 0.36/m~^. With bag constant 
ISOMeV the mixed phase region starts at density 0.13/m^^ and ends at 0.3/m~^. In figOl 
have separately plotted the EOS for Bg = ISOMeV showing the mixed phase region clearly. 
As it would be shown later that with only such choice of quark matter parameters I can 
attain the mass limit set by PSR J1614-2230. 

9 







1 1 1 




_ B =150 MeV 






g 






__. B^=160MeV 








1 , 1 





FIG. 6. Mass-radius curve with constant and varying bag pressure, having Bg = 170MeV. 
Assuming the star to be stationary and spherical, the Tolman-Oppenheimer-Volkoff 



(TOV) equations |32| gives the solution for the pressure P and the enclosed mass m. 



dP{r) _ G'm(r)e(r) [1 + P(r)/e(r)] [l + 47rr^P(r)/m(r)] 
dr l — 2Gm{r)/r ' 

^ = A.rh{r) , (12) 
dr 

G being the gravitational constant. Starting with a fixed central energy density e(r = 0) = 
ec, I integrate radially outwards until the pressure on the surface equals the one corresponding 
to the density of iron. This gives the star's radius R having gravitational mass 

pR 

Mg = m{R) =Att dr r^e{r) . (13) 
Jo 

For the NS crust, in the medium density range we add the hadronic EOS by Negele and 
Vautherin Jssj], and for the outer crust we add the EOS by Feynman- Metropolis- Teller 3^ 
and Baym-Pethick-Sutherland jssl . 

Fig [6] shows the gravitational mass M (in units of solar mass Mq) as a function of radius 
R, for constant and varying bag pressure Bg = 170 MeV. As the bag pressure varies and 
decreases towards the center of the star (at higher densities) the curve becomes stiffer as 
the effective matter pressure increases (bag pressure being negative). I find that a flatter 
EOS corresponds to a flatter mass-radius curve, and therefore the maximum mass of the 
star with varying bag pressure is higher than the non varying one. With such varying bag 
constant I plot the mass-radius curve with Bg = 160 MeV and 150 MeV (fig [7]). With the 
same qualitative aspect I find that the maximum mass of a mixed hybrid star obtained with 
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FIG. 7. Mass-radius curve with varying bag pressures, Bg = 160 MeV and 150 MeV. 



leOMeV is 1.84M0. The maximum mass with B„ = ISOMeV and m, = 300MeV, is 



2.01 solar mass. 



The discovery of high-mass pulsar PSR J1614-2230 [7] with mass of about 1.97M0, has 
set a stringent condition on the EOSs describing the interior of a compact star. They [3] 
quote the typical values of the central density of J1614-2230, for the allowed EOSs in the 
range 2no - 5no, whereas consideration of the EOS independent analysis of 361] sets the 
upper central density limit at IOhq. The maximum mass of a mixed phase EOS star with 
rris = ISOMeV is calculated to be 1.84 solar mass. The maximum mass for the mixed hybrid 
star can be increased to 2.01 solar mass, with = 300MeV having a varying bag pressure 
of Bg = ISOMeV. Only such choice of the quark matter parametrization can give rise to star 
which would satisfy the mass set by PSR J1614-2230. But with such choice of parameters 
the mixed phase region is small. This maximum mass limit is for this hadronic and quark 
matter EOSs. Stiffer EOS sets (like hadronic NL3 and quark ^ quark NJL model) for the 
mixed hybrid star can produce much higher maximum mass [37]. From the figure it is also 
clear that the maximum mass of the star corresponds to a radius of about 10km. Previous 
calculations have shown the maximum mass of a NS have radius greater than 12km, whereas 
the maximum mass of a SS corresponds to a radius of less than 9km. Therefore it is clear 
from my calculation that the mixed hybrid star has radius corresponding to the maximum 
mass, quite different from the neutron and strange star. They are not as compact as strange 
stars and their radius lies between the nuclear and strange star. 



11 



SUMMARY AND CONCLUSION 



In this work I have studied the maximum mass of a hybrid star having a mixed phase 
region. With the hadronic matter EOS having hyperons, and remaining in the simple 
MIT bag model I wanted to study what parameters value could give such high masses 
for a HS having a mixed phase region. The star has a dense quark core, a mixed phase 
intermediate region and hadronic outer region. The hadronic and quark matter EOS is 
simultaneously constructed according to relativistic mean field approach and MIT bag model. 
The mixed phase is determined in accordance with the Glendenning construction. All the 
phases are at chemical and mechanical equilibrium, and also they are charge neutral as a 
whole. With constant bag pressure Bg of 170 and ISOMeV (and rris = ISOMeV) I get EOS 
with considerable mixed phase region but with such parametrization the maximum mass 
of the star is about 1.5 solar mass. I therefore consider a density dependent bag pressure 
Bg, parametrized according to the Gaussian parametrization. The asymptotic value of the 
bag constant at high density is fixed at ISOMeV, which is its lowest value known from the 



experiments 



30l ]. With such varying bag pressure I can have a mixed phase region with 



Bg = 160MeV, but still the mass of the star is below 1.9 solar mass. To reach the mass 
limit set by PSR J1614-2230, for a mixed phase HS, I build the EOS with bag pressure of 
Bg = ISOMeV, having s-quark mass rris = 300MeV. For such choice of parameters values, 
the mixed phase region is small. Further lowering of bag pressure is not possible, as then 
the mixed phase disappears. The maximum mass for a mixed hybrid star with the given 
set of parameters is 2.01Mq. Another important results of my calculation is that the HS, 
with mixed phase, has radius (for the maximum mass) quite different from the neutron or 
strange star, their radius lying in between the neutron and strange star. 

After the discovery of PSR J1614-2230, setting the mass limit to 2 solar mass, new EOSs 



model has been proposed. Weissenborn et al. [381] showed that absolutely strange star can 
have mass above 2 solar mass is the effect of strong coupling constant and color supercon- 
ductivity is taken into account. Bednarek et al. 8| argued that EOS with hyperons having 
quartic terms involving hidden strangeness vector meson can reach such limit. Matsuda et 



al. 39j extended their calculation to hybrid stars, having a smooth crossover from hadronic 
to quark matter. For the mass to reach the maximum mass limit they showed that the 
crossover has to take place at low density and the quark matter has to be strongly interact- 
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ing. Using very stiff EOS sets (hadronic NL3 and quarlc quarlc NJL model) the maximum 
mass limit for the hybrid star can be raised much higher as shown by Lenzi & Lugones js^]- 
In my work, I also have shown that the maximum mass limit can be reached by a HS with 
mixed phase even with simple hyperonic nuclear matter EOS and MIT bag model quark 
matter EOS if I assume a relatively low density dependent bag pressure. 

Observationally the NS is characterised only by signals coming to us from its surface. 
Developments has been made on them to measure accurately the mass of compact stars but 
same cannot be done for their radius. Reasonable measurement of the radius of a compact 
stars could differentiate NS, SS and HS, as we have seen here that different EOS of matter 
gives different mass-radius relationship. As it is clear from my calculation and also from 
previous calculations that by suitable tuning of the parameters or by invoking new terms in 
the EOSs calculations the mass limit set by PSR J1614-2230 can be reached. Therefore to 
have a full understanding of the matter at extreme densities we need results not only from 
astrophysical observations but also from earth based experiments. 
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